Charge transfer between a superconductor and a hopping insulator 
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We develop a theory of the low-temperature charge transfer between a superconductor and a 
hopping insulator. We show that the charge transfer is governed by the coherent two-electron - 
Cooper pair conversion process, time reversal reflection, where electrons tunnel into superconductor 
from the localized states in the hopping insulator located near the interface, and calculate the 
corresponding interface resistance. This process is an analog to conventional Andreev reflection 
process. We show that the time reversal interface resistance is accessible experimentally, and that 
in mesoscopic structures it can exceed the bulk hopping resistance. 
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The transmission of the charge through the normal 
metal-superconductor interface occurs via the electron- 
hole conversion known as the Andreev reflection process: 
an electron incident from the metal side with an energy 
smaller than the energy gap in the superconductor is con- 
verted into a hole which moves backward with respect to 
the electron. The missing charge 2e (an electron has 
charge — e and a hole +e) propagates as an electron pair 
into the superconductor and joins the Cooper pair con- 
densate . Correspondingly, a Cooper pair transfer from 
the superconductor is described as the Andreev reflection 
of a hole. This Andreev transport channel is character- 
ized by the so-called Andreev interface contact resistance. 
Since transport current is introduced into a superconduc- 
tor via normal leads, the Andreev reflection phenomenon 
is a foundation for most applications of superconductors 
(see Ref. for a review). 

There exists however an important experimental sit- 
uation of the hopping insulator coupled to a measur- 
ing circuit via superconducting leads (see, for exam- 
ple, Q), where the conventional Andreev reflection pic- 
ture does not apply. The transport in hopping semicon- 
ductors occurs via localized [non-propagating) single par- 
ticle states with undefined momentum and therefore 
a Cooper pair on the superconductor side cannot form. 
A single particle transport through the interface is expo- 
nentially suppressed, cx e _A / T , where A is the supercon- 
ductor gap, the temperature, T, being measured in en- 
ergy units; therefore to explain the finite conductivity ob- 
served in experiments one needs Andreev-type processes 
capable to facilitate two particle transfer through the 
hopping insulator/superconductor interface allowing for 
Cooper pair formation. The possibility of such a transfer 
through the hopping-superconductor interface was dis- 
cussed in Ref. |5| but no quantitative theory of hopping 
transport - supercurrent conversion was presented. 

In this Letter we develop a theory for the transport 
through the hopping insulator-superconductor interface 



and derive the corresponding contact resistance. We 
show that the low-temperature charge transfer occurs via 
the correlated processes mediated by the pairs of hopping 
centers located near the interface. We demonstrate that 
this process resembles the conventional Andreev electron- 
to-hole reflection into a normal metal, the exponential 
suppression of transport specific to a single-particle pro- 
cesses being lifted. Thus, despite the limitation in the 
number of coherent hopping centers that can carry An- 
dreev transport, the resulting contact resistance can be- 
come low as compared to the resistance of the hopping 
insulator. However in mesoscopic structures the interface 
resistance can be comparable or even exceed the hop- 
ping resistance. The proposed mechanism resembles the 
so-called crossed Andreev charge transfer discussed 
recentlyin connection with a superconductor-dot entan- 
gler 0,0 ■ The difference is that in [1,0, H[ the transport 
mediated by artificial quantum dots was considered. In 
our case, the transport occurs via randomly located sites 
in the hopping insulator (HI), and the main problem one 
has to solve is finding the optimal configuration of the 
sites responsible for the charge transfer. Hereafter we 
will refer to the proposed charge transfer mechanism as 
to the time reversal reflection. 

Let a superconductor (S) and an HI to occupy the adja- 
cent 3D semi-spaces separated by a tunneling barrier (B). 
The presence of the barrier simplifies calculations which 
will be made in the lowest non-vanishing approximation 
in the tunneling amplitude Tq. This models the Schottky 
barrier usually presenting at a semiconductor-metal in- 
terface. In the linear response theory the conductance is 
determined by the Kubo formula 0] for the susceptibility, 



xH = i / ( i + (t),i(o) 



e lu3t dt 
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as Q = lim^^o u> 'Im^w). Here the current operator 



I(t) is defined as flpj ]: 

i(i)=iedT J d 2 r[a + (r,t)b(r,t) -h.c] , 

where r is the coordinate in the interface plane, a + (r,i) 
and b(r, t) are creation and annihilation operators in the 
semiconductor and superconductor, respectively, d is the 
electron localization length under barrier. The suscep- 
tibility, xi^i is calculated by analytical continuation of 
the Matsubara susceptibility 0| , 

X u(n)= [ (T T I(T)I(0))e iQT dr. (2) 
Jo 

Here T T means ordering in the imaginary time, (3 = 1/T. 
In the expression for (T t I(t)I(0)) one should expand to 
the second order with respect to the tunneling Hamilto- 
nian, 

H t (t) = dT J d 2 r [a+(r, r)6(r, r) + h.c] . (3) 

Keeping only those second order terms that contain 
(T T 6(r,r)6(r ,0)) (T T 6+(ri,ri)6 + (r 2) T 2 )) products and 
thus represent the time reversal scattering which we are 
interested in, one arrives at the expression 

(t t /(t)/(0)) = e 2 |T | 4 J dnd^Hd^^A + B): 

i 

A{i x i\) = F ( x - x o)F + (xi - x 2 )G(xi,x)G(x 2 ,x ) , 
B({xi}) = F(x - xi)F + (x - x 2 )[G(x ,x)G(x 2 ,xi) 
-G(x ,x 1 )G(x 2 ,x)], (4) 

where x = {r, r}, x = {r , 0}, x t = {n, n}; F(x - x') = 
(T T b(r, r)6(r', t')) is the anomalous Green function in the 
superconductor while G(x, x 1 ) = — (T T a(r, r)a + (r', t')) is 
the Green function in the hopping insulator. One can 
show that the Andreev-type process we are interested in 
is given by the first term of B({xi}) in Eq. (0J. The 
relevant diagram is shown in Fig. ^ Keeping only this 
term and using the Matsubara frequency representation 
one obtains 

XM (0) = 2Te 2 \T Q \ i d 4 f J] d 2 n ^ F(v - n, w n ) 

i ^ n 

xF+(r - r 2 ,u! n )G(r ,r,u! n - Q m )G(r 2 , r l5 -u) n )fi) 

where VL m = 27rmT and lo„ = (2n + l)wT. The 
normal Green's functions can be expressed through 
the wave functions of the localized states, <p s (r) = 
(tto 3 ) -1 / 2 exp(— |r — r s |/a), as 

G(ry,. n )=J2^^1. (6) 

s 

We have assumed that for all of the sites under consider- 
ation the voltage drops between the site and the super- 
conductor are the same. This is true when the partial 
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FIG. 1: The diagram describing the time reversal reflection. 
Lines with one arrow correspond to the Green's functions in 
the hopping insulator. They are associated either with the 
center 1, or with the center 2. Lines with two arrows corre- 
spond to anomalous Green's functions, see fill . Squares cor- 
respond to matrix elements of the tunneling Hamiltonian . 

interface resistance due to a time reversal pair is much 
larger than the typical resistance of the bond forming the 
percolation cluster. This situation resembles that consid- 
ered by Larkin and Shklovskii for the tunnel resistance 
between the hopping conductors |l2j . 

The anomalous Green function, F(R,u) n ), is 

F(R uj ) - [ - Jp-R/K 

7T 9m A sin(i?fc F ) _ W*»+»a 

= — e s . 7 

2 V 'A 2 + w 2 Rk F 

Here £ p = (p 2 — p 2 F )/2m, g m = mpF/^ 2 fi 3 is the density 
of states in a metal, kp = pf/^, while £ is the coher- 
ence length in a superconductor. Since F(R) oscillates 
with the period 27r/fc F integration over spatial coordi- 
nates along the interface yields the factor a 4 /k F \pi s \ 2 . 
Here pi s is projection of the vector R; s connecting the 
centers on the interface plane. Note that the dependen- 
cies on R and £ are similar to that given by Eq. (21) of 
the paper Q for the pair of the quantum dots near the 
superconducting interface. However the latter equation 
does not specify the dependences of the transmission co- 
efficients on the real physical parameters of the interface 
and the localized centers. 

The summation over the Matsubara frequencies, u n , is 
standard, 

T£/K) = /^/( £ )tanh^. 

The contour of integration closes the cuts |e| > A along 
the real axis. Upon analytical continuation one arrives 
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at the following expression for the conductance: 



ne*gl\T \W 
2hTk%a 2 



E 



7l(£,)n( £ ;)A 2 

|^| 2 (A 2 -^)' 



,-2(z s +zi)/a 



-2|pi. 



(8) 



Here n(e) = (e e ' + l) is the Fermi distribution, U c is 
the energy of the intcr-sitc Coulomb repulsion, and z-axis 
is perpendicular to the interface. 

In what follows we will replace by 
g 2 J d 3 ri d 3 r s dei de s where g is the the effective 
density of states in the hopping insulator. This is the 
density of states in the layer adjacent to the interface. 
Due to screening by the superconductor it is not affected 
by the Coulomb gap and can be considered as constant. 
Since we are dealing with the pairs close to the interface, 
the Coulomb repulsion is suppressed by screening. This 
screening can be conveniently regarded as an interaction 
of the charged particle with its image having the opposite 
charge. Thus the Coulomb correlation manifest them- 
selves as the dipole-dipole interaction and for p s i a 
one arrives at U c = e 2 a 2 / np%- Requiring it to me smaller 



than T one obtains a cut-off p s i > px = a(e / kclT) 



. , 1/3 . 

As a crude estimate, we take d A ~ fc^ 4 , while To 

„-A 



Tj,e~ A with T p ~ Bp. Bearing this in mind one finds 
9mTp/k F ~ 9m £ F/kp ~ 1- Since the ratio T p /(ak F ) 2 is 
of the order of the typical energy of the localized state, 
£d ~ h 2 /ma 2 , one arrives at the estimate 



G 9 — 2A 

gap T e d e 



Q n ~ -gaSe d e- 2A . 
ft 



(9) 



Here Q n is the conductance of a boundary between a nor- 
mal metal and a hopping insulator, while S is the contact 
area. The product gaSed is nothing but the number of 
localized centers within the layer of a thickness a near 
the interface and the factor gap T Sd expresses the prob- 
ability of finding a critical pair, i. e., a pair of nearly 
located hopping centers that dominates the time reversal 
reflection processes discussed above. 

The above approach holds, as we have already men- 
tioned, only if the resistance of the typical time reversal 
resistor (TRR) is much larger then that of the critical 
hopping resistor, i?/, = (h/e 2 -f) , where 7 is a dimcn- 
sionlcss factor depending on the mechanism of electron- 
phonon interaction and C is the hopping exponent |l3j . 
i. e., with the exponential accuracy, as long as 4A > (. 

There are many realistic situations where the barrier 
strength, A, is not too large; the Schottky barrier at the 
natural interface [l| is certainly the case like that. Conse- 
quently, if £ 3> 1, i. e., if the system is far from the metal- 
to-insulator transition point, the procedure of summation 
over the localized states should be modified. Namely, the 
choice of the pairs facilitating the charge transfer depends 
on the structure of the bonds connecting critical pairs to 
the rest of percolation cluster. 



According to the above considerations the voltage 
drops mainly on the bond connecting percolation cluster 
in the HI the critical TRR for which the distance between 
its pair components is less than the correlation length, C, 
of the backbone cluster. The incoherent electron trans- 
port can be ensured by a single one bond connecting the 
cluster to any of the TRR sites. Thus, the ratio Q/Q n is 
the probability to find a TRR contacting the percolation 
cluster. 

To estimate this probability let us consider the layer 
with the thickness of the typical hopping distance, 77,, 
near the interface where all the bonds of the backbone 
cluster have necessarily have a site within this layer. 
The total number of states in this layer is gSrhSh where 
Sh = T( is the width of the hopping band. This product 
can be estimated as ((3/8)(S/rh) 2 , where (3 is a numer- 
ical constant [l^. For the case of Mott variable range 
hopping (VRH), p « 20. 

The number of TRRs in this layer can be estimated 
as follows. Let us note first that the conserving energy, 
e s +Ei + U c , from ^-function in Eq. (JSJ) is associated with 
the band given by the broadening v = isq cxp(— 2rd/a). 
Here ra is distance to the nearest neighbor in HI. In- 
deed, the most natural source for the broadening of the 
resonance is coupling of the localized states. Secondly, 
since both electrons escape from the TRR through a 
single bond, the in-plane distance, p s i should not ex- 
ceed the typical distance between the hopping sites, 
rh = <C £. Keeping the exponential accuracy we ar- 
rive at the following criterion that the resistance of TRR 
is less than the resistance of a typical hopping resistor: 

4A+ln ^ + m ^N) + ^ + 2 (^ + ^) <c , (10) 

i/q T a a 

One may consider this equation as a generalization of 
the "connectivity criterion" to include the TRR. Here 
we deal with the independent variables e s , ej, z s and 
zi over which the averaging procedure should be done 
with an account of the restriction of Eq. I|10(l . Thus the 
number of the relevant TRRs is 



iTr 2 g 3 S I dei I de s / de d / r 2 dr d \ dz s dz/ / pdp 



x9 



max(e s ,e ; ) 2r d 2(z s + z t ) 



T 



(11) 



where a = 1 — [4A + h\{T / vq)] / Q < 1. Let us now mea- 
sure the energies in units of ash and lengths in units of 
ot/j, where Sh = TC,- Again, the product gr\eh can be 
estimated as /3/8, and we obtain the number of effective 
TRRs as Na ~ Aa 7 S/r 2 hl where 

A = 4ir 2 (l3/8) 2 f de, / de s / de d / if d dr d [ di h f d Vl 



x6 [max(e s ,e ; ) + i] d + r) s + m) - 1] « 0.1 



(12) 
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Consequently, Q/Q n ~ Aa 7 <C 1- One concludes that 
the difference between the "contact" resistances in nor- 
mal and superconducting states is dominated by the con- 
tribution of TRRs. Since Q N « R^S/C 2 , 

SR = g- 1 - g- 1 « -g- 1 = -R h (C 2 /SAa 7 ) . (13) 

Note that the interface resistance is of the order of the 
resistance of HI layer with the thickness ~ £/ (Aa 7 ) >> 
C. Since C ~ a£ 2 , one concludes that for £ > 10 the 
interface resistance can be comparable or even exceed 
the hopping resistance if the thickness of the sample (or 
of the contact) is < 10 fim. 

The resistance estimated above can be experimen- 
tally measured as a magnetoresistance in magnetic fields 
higher than the critical field for superconductivity (sim- 
ilar effect for quasiparticle channel was studied in |5j). 

We were implicitly assuming so far that the variable 
range hopping occurs according to the Mott's law. This 
assumption certainly holds near the interface where the 
Coulomb gap is screened by a superconductor. However 
in the bulk of HI the Efros-Shklovskii (ES) law (3 
can become the dominant hopping mechanism. Then the 
value of C in the connectivity criterion l|10f) will be con- 
trolled by the Coulomb gap, C — > Cer = {Pie 2 / naT) 1 / 2 , 
where (3\ is a numerical constant [ujQ]. In this case, 
it turns out that each bond of the ES backbone cluster 
finds some TRR ensuring charge transfer. Thus in the 
limiting case of a weak tunneling interface barrier the 
contact resistance will be the same for both normal metal 
and superconductor leads. This fact can be used to dis- 
criminate between Mott and Efros-Shklovskii laws in the 
situation when it is difficult to do so from temperature 
dependence. 

Note that, in principle, the charge transfer involving 
double occupied localized states is possible. However, 
such a process would require an additional activation ex- 
ponential factor, cx e~ u ^ T , where U is the on-site corre- 
lation energy. One can also consider processes where a 
double occupied center (so-called D~-center) serves as an 
intermediate state for the phonon-assisted two-electron 
tunneling. This channel is unfavorable (at least in the 
case of a large interface barrier) because of the above- 
mentioned exponential factor and a small pre-exponential 
factor due to phonon-assisted tunneling. For the weak 
tunnel barriers the conductance is controlled by "typical" 
hopping sites. In this case D~ channel is suppressed ci- 
ther by the additional tunneling exponential, cx e _4r,1 / a , 
or by a small probability to form a close triad of hopping 
sites. Therefore the D~ channel can be also neglected. 

To conclude, we have developed a theory of the low- 
temperature charge transfer between a superconductor 
and a hopping insulator and calculated the interface re- 
sistance. This resistance is dominated by time reversal 
reflection processes involving localized states in the insu- 
lator. It is the time reversal reflection process that allows 
the low-temperature measurements of hopping transport 



utilizing superconducting electrodes in the experimental 
setups. In the Efros-Shklovskii VRH regime, the corre- 
sponding interface resistance is small as compared to the 
bulk hopping resistance and is nearly equal to the resis- 
tance at the interface between the HI and normal metal. 
On the contrary, in the Mott hopping regime (relevant, 
in particular, for 2D gated structures), the interface re- 
sistance grows much larger and becomes commensurate 
(or even exceeds) to the bulk hopping resistance. This 
effect is especially pronounced in the mesoscopic sam- 
ples. The contribution from the interface resistance can 
be detected by application of the external magnetic field: 
the relatively weak magnetic field will drive the super- 
conductor into the normal state, but will not affect the 
hopping transport eliminating thus time reversal reflec- 
tion process. This effect holds even in the case where the 
interface contribution is less than the typical resistance 
of the hopping system itself. 
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